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1. Introduction 

Quenched randomness may have a profound effect on the nature of phase transitions. If 
there is a continuous phase transition in the perfect system then, according to the Harris 
criterion 0, the relevance of the perturbation is connected to the sign of the specific heat 
exponent a in the pure system. The two-dimensional random Ising model with a = 
represents the borderline case of the perturbational theory. Indeed, that model has been 
the subject of intense investigations to clarify its critical properties 0, [|. 

According to field-theoretical studies [0, || the randomness is, in the renormalization 
group sense, a marginally irrelevant perturbation, therefore it leads to logarithmic 
corrections to the power-law singularities of the pure model. For example, the bulk 
magnetization, m^, and the correlation length, £, are expected to behave near the 
transition point as 



t^llntr 1 / 16 , (1) 



and 



e-r^lntl 1 / 2 , (2) 



where t = \T C — T\/T c is the reduced temperature. The critical spin-spin correlation 
function G(r) averaged over several samples has a pure power law decay || 

G(r) ~r~ 1/4 [A + B/(\nr)' 2 ] , (3) 

whereas the typical correlation function calculated in a large single sample is 
conjectured^ to decay as 

G(r) ~r" 1/4 (lnr)" 1/8 . (4) 

The above conjectures are found to be in agreement with numerical results of large- 
scale Monte Carlo (MC) simulations ]7|, (|, || and transfer matrix calculations [[|, [10], [TT 



However, also conflicting interpretations of the numerical findings have been suggested, 
invoking dilution-dependent critical exponents and weak universality |T3, 



In this paper, we consider the boundary critical behaviour of the two-dimensional 
random bond Ising model. The surface critical properties of the perfect model are 



exactly known since many years [pL4| . For example, the asymptotic behaviour of the 
surface magnetization, m 1; and the correlation length, £y, measured parallel to the 
surface, is given by 

mi ~ t 1/2 (5) 

and 



£\\~t 1 , 



(6) 
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whereas the critical surface spin-spin correlation function has the asymptotic decay: 

G s (r) ~ r" 1 . (7) 

Thus the corresponding critical exponents are f3\ = 1/2, u» = 1 and rj\\ = 1. No field- 
theoretical results are available for the random case. However, it seems reasonable to 
expect, in analogy to the bulk properties, that the randomness is a marginally irrelevant 
variable at the surface fixed point as well. Then one might obtain logarithmic corrections 
to the asymptotic behaviour of the perfect model. 

In the present study, we performed extensive numerical investigations to illuminate 
this issue by determining the surface critical properties of the Ising model with nearest- 
neighbour random couplings on the square lattice. In our first approach, we used 
large-scale MC techniques and computed the surface magnetization and the complete 
magnetization profile of the model. Our second method is based on the star-triangle 
(ST) transformation. By that method we calculated both the surface magnetization 
and the surface correlation function of the model. By the two, in several respects 
complementary approaches, we determined numerically the complete set of surface 
critical exponents, including the surface magnetization exponent /3i, the correlation 
length exponent v\\ and the decay exponent of the critical suface correlations m. Note 
that some of the MC results on the surface magnetization have already been announced 



in a short communication 15 



The paper is organised as follows. The MC results on the surface magnetization 
and the magnetization profiles are presented in Section 2. In Section 3 we describe the 
ST approach as applied to the random Ising model and discuss the numerical results on 
the surface magnetization and the surface correlation function. The main conclusions 
are given in the Summary. Some details of the ST method have been transferred to the 
Appendix. 



2. Monte Carlo simulations 

Let us consider the Ising model with nearest-neighbour ferromagnetic couplings, where 
the spins Sjj(=±l) are situated on the sites of a square lattice. A surface may be 
introduced by cutting the coupling bonds along one of the axes of the lattice, leading 
to the (10) surface, or along the diagonal, leading to the (11) surface. In the MC 
simulations, we studied systems with two parallel surface lines, each line having L sites. 
Each row perpendicular to the surface consists of K sites. The spins in the first and last 
row are assumed to be connected by periodic boundary conditions. The lines parallel to 
the surfaces are numbered by the index i, i.e. % = 1 and i = K denote the two surface 
lines. The index j refers to the position along a line, running from 1 to L. The total 
number of spins is K x L. The aim of the simulations is to determine thermal properties 
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of the semi-infinite system, where K, L — ► oo; therefore finite-size effects need to be 
studied with care. 

The interaction between neighbouring spins may be either 'strong', J\ > 0, or 'weak', 
< J2 < J 1 ■ Strong and weak couplings are distributed randomly, with p (or f — p) 
being the probability of a weak (or strong) bond. If both interactions occur with the 
same probability, p = 1/2, then the model is self-dual ]TJ|. The self-dual point is 
located at 

tanh( JjksTc) =exp(-2J 2 /k B T c ). (8) 

determing the critical temperature, T c , of the bulk Ising system (K, L — > 00, and full 
periodic boundary conditions), if the model undergoes one phase transition. Indeed, 
results of previous simulations support that assumption [|, [7]. The simulations were 
done for the self-dual case, i.e. at p= 1/2. 

Certainly, one expects that both bulk and surface will still order at the bulk critical 
temperature, T c , in a two-dimensional model with short-range interactions: The one- 
dimensional surface does not support any separate ordering, so that one encounters the 



'ordinary transition' [17], |18|| . 

Varying, in the self-dual case, the ratio of the strong and weak couplings, r = J%j Ji, 
one may change the degree of dilution. At r = 1, one recovers the perfect Ising model, 
while r = corresponds to the percolation limit, where T c = 0. As had been shown 
before, the crossover to the randomness dominated bulk critical behaviour may be 
monitored conveniently by choosing r in the range of 1/10 to 1/4. Then the crossover 
length, at criticality, ranges from a few to about 20 lattice spacings |7[]. Indeed, we 
simulated the random model at these two values, r = 1/4 and 1/10, augmented by 
computations for the perfect model, r = 1. 

Most of the simulations were performed for the (11) surface, albeit a few runs were 
also done for the (10) surface to compare with exact results. For the (11) surfaces, we 
usually set L = K/2, with K ranging from 40 to 1280 to check for finite-size effects. For 
the (10) surfaces, quadratic systems were studied. We averaged over an ensemble of bond 
configurations (or realizations). The number of realizations typically ranged from at 
least 15 for the largest systems up to several hundreds for the small systems. In general, 
the one-cluster flip Monte Carlo algorithm was used (mainly for testing purposes, 
we also applied the single-spin flip method), generating, close to the critical point, 
several 10 4 clusters per realization. Note that the statistical errors for each realization 
were significantly smaller than those resulting from the ensemble averaging. To avoid 
inaccuracies due a, possibly, unfortunate choice of the random number generators, we 
compared results obtained from shift register and linear congruential generators. 

The crucial quantity, computed in the MC simulations, is the magnetization profile. 
It is described by the magnetization per line, m{i) =< \ s i,j\ > / L, where s^j denotes 
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the spin in line i and row j, with % — 1,2, ...K, and summing over j = 1,2, ...L. The 
absolute values are taken to obtain a non-vanishing profile for finite systems, as usual. 
The surface magnetization is given by mi = m(l) = m(K). 

Because the distribution of the random bonds is the same in the bulk and at the 
surface, one may expect a monotonic decrease of m(i) on approach to the surface, due 
to the reduced coordination number at the surface (being two for the (11) and three 
for the (10) surface). This behaviour is illustrated in Fig. 1, comparing magnetization 
profiles of the perfect, r = 1, and random, r = 1/4, Ising model with a (11) surface, 
at the same distances from T c , measured by the reduced temperature t — \T — T c \/T c . 
The critical point, T c , follows from (1). Obviously, randomness tends to suppress the 
magnetization, at fixed value of t. The profiles display a pronounced plateau around the 
center of the systems, at which the bulk magnetization, 7m,, is reached. The existence of 
the broad plateau indicates that the linear dimension K of the MC system is sufficiently 
large to compute, for instance, the surface magnetization of the semi-infinite system. Of 
course, in addition one has to monitor possible changes of mi with L, to have possible 
finite-size effects due to that dimension under control.- Note that is known exactly 
in the perfect case [fL4| , and very accurately in the random case 0. 



For the perfect two-dimensional Ising model with a (10) surface, the complete 
magnetization profile has been calculated exactly in the continuum limit ]T^, [2tJ. 
In particular, the profile approaches the bulk value in an exponential form, with 
rrib — m(i) oc exp(—i/^ r ), where £ r is the bulk 'correlation range', which becomes 
only asymptotically, as T — > T c , identical to the bulk 'true correlation length' ET 



Indeed, we tested the accuracy of our simulational data by comparing them, for the 
(10) surface, to the exact expression. In addition, we found that the same correlation 
range determines the exponential approach of the magnetization towards its bulk value 
in the (11) case as well. 

For the perfect two-dimensional Ising model with a (11) surface, exact results exist 
for the surface magnetization, mi, and the magnetization in the next line, m(2) ||22|| . 
Again, the Monte Carlo data, obtained with modest computational efforts, agreed very 
well with the exact results, as shown in Fig. 2. In the figure, the 'effective exponent' 
(3{i) e ff is depicted, defined by 

fi(i) eff (t) =dln{m(i))/dln(t) (9) 

Certainly, as t — > 0, the effective exponent acquires the true asymptotic value of 
the critical exponent f3(i). For example,the asymptotic critical exponent of the surface 
magnetization is f3{l) = Pi = 1/2, being, by the way, identical for (11) and (10) surfaces. 
Because the magnetization m(i) is computed at discrete temperatures t k , we use in 
analysing the simulational data, instead of (9) the corresponding difference expression, 
with t = {t k + t k+ i)/2. The error bars, included in Fig. 2, have been calculated in a 
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conservative fashion, getting the bounds for (3 by comparing the upper (lower) limit of 
m(i) to the lower (upper) limit of m(i + 1), where the bounds of the magnetization 
are computed in the standard way from the ensemble averaging. Alternately, we 
also computed the error bars from usual error propagation, which turned out to be 
appreciably smaller. 

In Fig. 2, the temperature dependence of the effective exponent fl(i) e ff deeper in the 
bulk is also displayed. For example at i = 10, one readily observes the crossover from the 
bulk effective exponent (as follows from the exact expression for the bulk magnetization 
|14| ) to the surface dominated behaviour, when the correlation length becomes large 
compared to the distance from the surface. In general, at finite and arbitrarily large 
distances to the surface, /3(i) e ff will always converge, on approach to T c , to the surface 
critical exponent, (3i = 1/2, and not to the bulk critical exponent, (3 =1/8. Analogous 



observations have been reported for three-dimensional Ising models with surfaces [23 



The main aim of the Monte Carlo study has been to estimate (3\ in the random case. 
Results of the extensive simulations are summarized in Fig. 3, depicting the effective 
exponent /3(l) e //(t) at r = 1/4 and r = 1/10, compared to its exactly known form for 
the perfect case, r = 1. Typical error bars, increasing closer to criticality, are displayed, 
based on standard error propagation resulting from the variance in ensemble averaging 
of mi(tfc) and mi(tk + i). Data obviously affected by finite-size effects have not been 
included in the figure. 

As seen from Fig. 3, at fixed distance from the critical point, t, /3(l) e // rises 
systematically with increasing dilution, reflecting the decrease in m x with stronger 
randomness. However, asymptotically, t — > 0, it is well conceivable that the surface 
critical exponent will coincide in the perfect and dilute cases, with (3\ = 1/2. Indeed, a 
reasonable estimate, both for r = 1/4 and r = 1/10, is (3\ =0.49 ±0.02. 

Thence, the simulations demonstrated that the critical exponent fli is rather robust 
against introducing randomness simultaneously in the bulk and at the surface. Note that 
f3\ remains 1/2 too, when only the surface bonds of the two-dimensional Ising model are 
randomized as described above, but keeping a unique bulk coupling, as we confirmed in 
simulations. Interestingly enough, in the three-dimensional case, introducing random 
nearest-neighbour strong and weak surface bonds, but having only one interaction for 
the bulk couplings, seems to be an irrelevant perturbation as well, i.e. the surface critical 



exponent seems to be the same as for the perfect surface, (3\ ~ 0.80 [23]. This robustness 



may indicate that the bulk critical fluctuations play a crucial role for the surface critical 
exponent, albeit it is not determined by bulk critical exponents [17 , |18 . If that is true, 
then our result for the two-dimensional case with random bulk and surface interactions 
suggests that the bulk critical fluctuations are not very sensitive towards dilution (in 
accordance with the theory of, at most, logarithmic modifications of the asymptotic 
power-laws describing critical behaviour of the perfect system in two-dimensional Ising 
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models [||, |J). We shall come back to this aspect in the next Section. 
3. Star-triangle transformation 

The star-triangle transformation was introduced by Hilhorst and van Leeuwen[^4j], and 



used later by others |25|, |27| to calculate the surface magnetization and the surface 
correlations in layered triangular lattice Ising models. Here we generalize the method 
for non-translationally invariant systems. 

3.1. Star-triangle approach to boundary behaviour 

The method is based on an exact mapping of the original triangular model, with 
couplings {Ki}, {K 2 } and {-^3}, to a hexagonal model with couplings {pi}, {P2} and 
{^3} denoted by dashed lines in Fig. 4. In the transformation the right-pointing triangles 
are replaced by stars such that the couplings are related by 

K 1= *] n ( cosh (Pi +P2+ Ps) cosh(-p! +p 2 + p 3 ) \ / 1q n 
1 4 \ cosh (pi +p 2 -ps) cosh (pi - p 2 + p 3 ) / 

and its cyclic permutation in the indices i = 1,2, 3. In the second step of the mapping 
the left-pointing stars of the hexagonal lattice are replaced by triangles resulting in 
a new triangular lattice, which is denoted by dotted lines in Fig. 4. Iterating this 
procedure a sequence of triangular Ising models is generated (n = 0,1,2,.. .) from the 
original model with n = 0. 

As seen in Fig. 4, the surface spins of the n-th and the (n + l)-th models are 
connected by the surface couplings of the intermediate hexagonal model. In this 
geometry, the thermal average of the Z-th surface spin of the n-th model, denoted by 
(s\ n ^) = o~i is connected to the thermal averages of the neighbouring spins and 
of the (n+l)-th model, where the corresponding surface couplings of the hexagonal 
lattice are denoted by and As shown in the Appendix, one has 

al n) = ak^a^ + at +1) o-£ +1) , (11) 



where 



1 - tanh 2 (p{ ra h 
l-tanh^Vanh 2 ^ * 



$ - tanhfa?) , (12) 



while in one should interchange pf^ and p^ . Now using the vector notation 
for the surface spin exceptational values in the n-th model and similarly a < - n+2 ' ) for the 
(n + 2)-th model we obtain the relation 



a 



(n) = A (n+l) A (n+2) a (n+2) ^ 13 ^ 



where the non- vanishing elements of the A^ n+1 ^ matrix are given by a+ +1 \l) and 
a^L +1 \l) in terms of the surface couplings of the (n + l)-th hexagonal lattice, equation 
(P"2|), and similar relation holds for A^ n+2 \ (We consider two successive steps in ( p"3D in 
order to avoid complications with the different parity of the odd and even number of 
transformations.) Now taking the boundary condition lim^oo = (1, , 1, . . . , 1) we 
obtain for the average surface magnetization 

1 L 

mi = lim - J2 <7i = iim f(n) , (14) 



with 



l-oo L l=1 



1 - 
fin) = lim — 

*J=1 



(15) 



We note that f(n) in (JT3J) is formally equivalent to the partition function of an n- 
step directed walk (polymer) in a random environment, where the (random) fugacities 
corresponding to the k-th step of the walk are contained in the A^ k ' matrix, which is 
just the transfer matrix of the directed walk. 

Next we consider the average connected surface correlation function defined as 

G s (l) = lim I]T[<5 i+ ^> - (s i+l )( Si )] . (16) 

i=l 

As shown in the Appendix, the surface correlations in the n-th triangular model are 
connected to those in the (n + l)-th model, and the relation is given in terms of the 
surface couplings of the intermediate hexagonal lattice, equation (30), similarly to (Tfl|). 
Furthermore, as we argue in the Appendix, in the asymptotic limit (/ 3> 1) the surface 
correlation function can be expressed by the partition function f(n) of the corresponding 
directed walk, 

r°° I ( I 2 \ 

G.(J)oc/_ ^-^exp(--)[/ 2 H-/ 2 (oo)]. (17) 



n 



Thus the surface properties of the model are connected to the asymptotic behaviour of 
f(n) in (p"5"D. For different temperatures, corresponding to different thermodynamical 
phases of the random Ising model, the f(n) function has different asymptotic behaviour, 
as can be seen in Fig. 5 for a dilution of r = 1/10. 

In the ordered phase, T < T c , f{n) approaches a finite limit, the surface 
magnetization m 1; through an exponential decay, 

f(n) =m 1 (T) + Aex P (-n/Zj) T < T c . (18) 

For T > T c the limiting value of f(n) is zero, which corresponds to vanishing surface 
magnetization, and the decay for T > T c is exponential, 

f(n) ex exp(-n/e 2 ) T > T c , (19) 
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whereas at the critical point, it has the form of a power law 

f(n) (xn-< T = T C . (20) 



We argue that £y in equations ( |T8| ) and (|T9| ) is the surface correlation length, below and 



above the critical point, respectively. Indeed, substituting (^) or ( |T9| ) into ( |TTD and 
evaluating the integral by the saddle-point method, we get 

G s (Z)cxexp(-Z/£||), (21) 

in accordance with the definition of the surface correlation length. 

At the critical point, where f(n) as in fl20|), the surface correlation function in 
equation ( |TTD leads to a power law decay G s (l) ~ Z~ 47 . Thus the decay exponent, r)\\, of 
the critical surface correlation function is given by 

77,1 = 4 7 . (22) 

We conclude at this point, that we have obtained a complete description about 
the surface properties of the random Ising model by the star-triangle method. In 
the following, we shall use the above formalism to calculate numerically the surface 
magnetization, the critical surface correlations, and the correlation length. 



3.2. Numerical results 

In the actual calculations, we considered the random Ising model of the MC simulations, 
with a (1,1) surface, on a strip of width L of a diagonal square lattice (which can be 
considered as a triangular lattice with vanishing vertical couplings) and imposed periodic 
boundary conditions. To reduce finite size effects, we considered relatively large strips 
(with L up to 512). We checked that the difference between the results for the two 
largest strips (L = 256 and L = 512) is essentially negligible, doing up to n = 2000 
iterations |28|]. We calculated the partition function f{n) as a function of n, averaging 
over several (typically around twenty) random configurations of the couplings. The ratio 
r between the two, weak and strong, random couplings was chosen to be 1, 1/4, and 
1/10, as in the simulations; both couplings occur with the same probability, p = 1/2. 

We start with the analysis of the results in the ordered phase, i.e. T < T c . For 
a given temperature, f{n) approaches the surface magnetization mi, see (|T4]), which 
is found to agree (within the error of the calculations) with MC data presented in 
the previous Section. Approaching the critical point, the convergency of f(n) with n 
becomes slower, in accordance with the form of the correction term in (|18D . Accordingly, 
accurate estimates become more difficult. As in the case of the simulations, the ensemble 
sampling over different configurations seems to be, however, the main source of error. 

From the values for mi(t) at different reduced temperatures t, we determined 
effective surface magnetization exponents /3(l) e ff(t), as defined in (|J). The estimates 
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of the effective exponents obtained from the star-triangle method are close to those 
found by the MC technique, see Fig. 3. Thus we confirm that f3i is rather robust against 
introducing randomness in the two-dimensional Ising model. 

At T c , we studied the surface correlation function, as follows from the partition 
function /(n). As shown in Fig. 6, f(n) exhibits, with n ranging from 100 to 1000, on 
a log- log plot (ln(/(n)) vs. ln(ra)), an almost linear behaviour. The average slope then 
defines an average decay exponent j av , see Q2"(]|). For the perfect model, our estimate 
agrees nicely with the exact value ^ pure = 1/4. In the random case, the average exponent 
decreases with rising randomness, i.e. decreasing value of r. For 100 < n < 1000, we 
obtain the average exponents ^ av = 0.228 and ^ av = 0.207, at r = 1/4 and r = 1/10, 
respectively. Based on these estimates, one may argue, that also the decay exponent 
97||, see equation (|22|) , varies with dilution, r. However, a more detailed analysis is 
needed to substantiate or repudiate these statements. For instance, looking at the local 
effective exponent, defined in analogy to @, a slight increase of that exponent with 
increasing n is observed. Indeed, the data for f{n) depicted in Fig. 6, may be fitted 
by the power law of the perfect model modified by logarithmic corrections with almost 
identical confidence (doing least square fits) as by the power laws with the average, 
dilution dependent exponents. Thence, our data leave room to different interpretations. 

In the disordered phase of the model, T > T c , we studied the correlation length 
from the asymptotic decay of f(n) in (|i~9|). Examples of the results of our calculations 
are shown in Fig. 7, plotting ln(/(n)) as a function of n at several temperatures t for 
r = 1/10. As seen from that figure, f(n) seems to exhibit an exponential decay, with 
£||(t)~ 2 following from the slopes of the curves. Approaching the critical temperature 
T c , the correlation length is expected to diverge as £y(t) ~ t~ v . From data at t > 0.05, 
we calculated average critical exponents z/ a „(t), with v av = 1.07(2) at r = 1/4, and 
v av = 1.13(6) at r = 1/10. These average exponents are larger than the asymptotic 
exponent of the pure model, v pur e = 1, see (||), and vary with the degree of dilution. Note 



that similar values have been obtained before by using finite size scaling |TTJ. However, 
those average values have been interpreted as reflecting logarithmic corrections to the 
power law of the perfect case [|lTJ. Again, we cannot rule out that possibility. 



For further interpretation of our data, we consider the scaling relation [17 



A = ^||/2, (23) 

which is satisfied, within the errors of the estimates, by the average exponents, both 
for r = 1/4 and r = 1/10. Following the alternate interpretation involving logarithmic 
corrections, the critical surface correlations, described by r)\\, would be then affected by 
logarithmic corrections, due to the correction terms in the correlation length (and their 
presumed absence in the surface magnetization). The above scaling law, (|23|), can be 
obtained by relating the surface correlation function between two spins at a distance of 
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the correlation length, £(£), to the square of the surface magnetization, 

G,(t(t)) ~ rn\{t) , (24) 
in the limit t — > (when logarithmic corrections are present, such a relation has been, 



for instance, checked for the q = 4 state Potts model p9|). Then, supposing logarithmic 



terms in the surface correlations, but not in the surface magnetization, one easily arrives 
at the conjecture 

G a (l) ~Z- 1 (lnZ) 1 / 2 . (25) 



4. Summary 

In this paper, the boundary critical properties of the two-dimensional random Ising 
model have been studied by MC techniques and by the star-triangle (ST) approach. In 
the simulations, we computed magnetization profiles, allowing to monitor surface and 
bulk quantities simultaneously On the other hand, by the ST method we calculated the 
surface magnetization as well as surface correlation functions. Both methods provide 
data for the surface magnetization which are in very good agreement, demonstrating 
the correctness and accuracy of the two approaches. 

To analyse the behaviour of the random Ising model in the critical region, we 
considered three singular quantities: the surface magnetization, the (surface) correlation 
length and the critical surface correlation function. The surface magnetization of the 
dilute model, as obtained from both methods, follows closely the power law of the 
corresponding perfect model, where (3\ = 1/2, showing the robustness of that exponent 
against even fairly strong randomness. 

The behaviour of the other two singular quantities, the critical surface correlation 
function and the correlation length, as determined from the ST method, is rather subtle. 
Note that in the ST method we averaged over the logarithm of the surface correlation 
function, leading to information about the typical behaviour of the correlation length. 
The numerical estimates for the critical exponents of both quantities, i.e. u\\ and 7711, are 
found to be dilution dependent, when calculating average exponents. Similar findings 
have been reported before for bulk exponents in the two-dimensional random Ising 
model. These non-universal average bulk exponents have been interpreted previously 
either as reflecting the true asymptotics (satisfying weak universality) or as being due to 
logarithmic corrections of the power laws in the perfect model (in accordance with field- 
theoretical predictions). Our numerical data for the surface quantities leave room to 
both types of interpretation, as concerns bulk and surface critical properties. Extending 
the field-theoretical predictions and attributing the apparent variation of the average 
exponents with the degree of dilution to logarithmic corrections, we conjectured, in 
equation fl2"o]), the asymptotic form of the critical surface correlation function. 
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Appendix 



We consider the first two layers of a hexagonal Ising lattice (Fig. 4), where a surface spin 
s is connected to the second layer spins s + and s_ by couplings p + and p_, respectively. 
We are interested in a relation between the thermal averages (s+), (s_) and (s). 
We start by considering the conditional probability 

exp(p + ss + +P-SS-) 



Pis) 



(26) 



s+ ' s ~ J2 S exp(p + ss + + p^ss^j 
with fixed values of s + and s_. Under this condition the expectational value of s is 
given by 

(8)1. . 



tanh(p + ss + + 



tanh 



+ s_ 



P+-P- 



(27) 



which can be evaluated using the fact that s. 

tanh(p + + p_) + tanh(p + — p_) 



±1 and s_ = ±1 as 



S+,S- 



tanh(p + + p_) + tanh(p_ — p^ 



+ 8. 

2 

At this point, one can sum over the variables s + and s_ leading to 
. tanh(p+ + p_) + tanh(p + - p) 

\ s ) = \ s +) o 



(2? 



+ (0 



tanh(p + + p_) + tanh(p_ — p + ) 



(29) 



which is equivalent to equation (|TT|). 

The connected surface correlation function of the n-th triangular model (i+l, i) 
(si+i s i) ~ ( s i+i)( s i ) an d that °f the (n + l)-th model are related by 

gW(i + /,<) = a^ + att +1) 9 (n+1 \i + / + l,i + 1) + a^+oi-4- +1 ^ ( " +1) + I + 1, < - 1) 

Z-l,i+l) 



(30) 

which can be obtained along the lines of (|TTD . Iterating the expression in (30), one 
obtains a sum, each term of which can be formally represented by two directed walks, 
which start at positions i + I and i, respectively. If the two walks meet at step n and 
at some position j, then g^ n \j,j) = 1 — (Sj) 2 and the walks annihilate each other. In 
the n — > oo limit, the non- vanishing contribution to g^°'(i + l,i) = g(i + l,i) is given 
by those processes, which are connected to annihilated walks. In the transfer matrix 
notation the average surface correlation function is given by 

i l 



Gs(l) 



L 



EE 

j=l a.w. 







n 






_m=l 


i+l, j 


_fc=l 







(31) 



13 



The asymptotic behavior of this expression can be obtained by noticing that the 
transverse fluctuations of directed walks have Gaussian nature, and the corresponding 
probability distribution is sharp. Consequently, in the large I (and large n) limit it is 
enough to consider the typical processes. Then there are two factors in (|31|), which are 
both approaching the partition function of n-step directed walks, f(n), see (fH)|), and 
these contributions should be multiplied by -P n (0> the ratio of those walks which are 
annihilated at the n-th step. In this way, we obtain 

G s (l)^E P n(mf 2 (n)-f(oo)], (32) 

n 

which in the continuum approximation is given in fll7l). 
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Figure captions 



Figure 1. Magnetization profiles m(i) of two-dimensional perfect (squares) and 
random, r = 1/4 (circles), Ising models with (11) surfaces, at t = 0.2 (open symbols) 
and t = 0.05 (full symbols). Systems of size 160x80 were simulated. 



Figure 2. Effective exponent /3(i) e //, with i=l,2,3, and 10, vs. reduced temperature t, 
for the perfect Ising model with (11) surface. The solid lines denote exact results[14,22]. 
Monte Carlo data for systems of sizes 80x40 (t > 0.3), 160x80 (0.07 < t < 0.3), and 
320x160 (t < 0.07) are shown. 



Figure 3. Effective exponent /3(l) e //, vs. reduced temperature t for the random two- 
dimensional Ising model with (11) surface, at r=l/4 (circles) and r=l/10 (triangles). 
Systems of sizes 80x40 (t > 0.3), 160x80 (t = 0.275), 320x160 (0.1 < t < 0.275), 
640x320 (0.05 < t < 0.1), and 1280x640 (t < 0.05) were simulated. The solid line 
denotes the exact result in the perfect case[22]. 



Figure 4. Mapping of the original triangular lattice (solid line) to an equivalent 
hexagonal lattice (dashed line) and further to a new triangular lattice (dotted line) 
using the star-triangle transformation. The surface spins of the n-th model (Si ) and 
those of the (n + l)-th model (S^ +1 \ S^" ' ) are connected by the surface couplings 
of the intermediate (n + l)-th hexagonal model (p;? , P;™ )■ The couplings {Ki}, 
{Pi}> i = 1,2,3 appearing in the star-triangle relation in ( |i~0| ) are also indicated. 
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Figure 5. Magnetization of the (1,1) surface of the two-dimensional random 
Ising model with r = 1/10 as a function of the temperature. The finite iteration 
approximants f(n) of the star-triangle method in equation ( [l5| ) are indicated by circles 
(n = 128), squares (n = 256), triangles (n = 512) and by crosses (n = 1024). The 
asymptotic behaviour of fin) is different for T < T c , T > T c and at T = T c , as given 
in (|l|),© and ©, respectively. 



Figure 6. Finite iteration approximants to the surface magnetization, f(n), as a 
function of n in a log- log plot, at the critical point of the two-dimensional random 
Ising model with dilution r = 1/10 (squares) and r — 1/4 (circles), compared with 
the perfect model (triangles). The slope of the curves, indicated by straight lines, is 
related to the average decay exponent 7711 of the critical surface correlations through 
@, see text. 



Figure 7. Finite iteration approximants to the surface magnetization, f(n), as a 
function of n in a semi- logarithmic plot, at different reduced temperatures t = 0.1 
(triangles), t = 0.2 (circles) and t = 0.3 (squares) above the critical point of the two- 
dimensional random Ising model with r = 1/10. The slope of the curves, indicated by 
straight lines, corresponds to the inverse square of the average correlation length, see 
equation jl9|). 



